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We report on the nonlocal gauge invariant operator of dimension two, Fpv(£'^) '/^v We are able to localize 
this operator by introducing a suitable set of (anti)commuting antisymmetric tensor fields. Starting from 
this, we succeed in constructing a local gauge invariant action containing a mass parameter, and we prove 
the renormalizability to all orders of perturbation theory of this action in the linear covariant gauges using 
the algebraic renormalization technique. We point out the existence of a nilpotent BRST symmetry. Despite 
the additional (anti)commuting tensor fields and coupling constants, we prove that our model in the limit of 
vanishing mass is equivalent with ordinary massless Yang-Mills theories by making use of an extra symmetry in 
the massless case. We also present explicit renormalization group functions at two loop order in the MS scheme. 
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I. INTRODUCTION 

We shall consider pure Euclidean SU (N) Yang-Mills theo- 
ries with action 



with 



(4) 



Sym — 



pa pa 



(1) 



where A^, a = 1, ...,A^ — 1 is the gauge boson field, with as- 
sociated field strength 



(2) 



The theory Q is invariant with respect to the local gauge 
transformations 



5a;; = co^ 



(3) 
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denoting the adjoint covariant derivative. 

As it is well known, the theory Q is asymptotically free 
Qll^, i.e. the coupling becomes smaller at higher energies 
and vice versa. At very high energies, the interaction is weak 
and the gluons can be considered as almost free particles. 
However, in spite of the progress in the last decades, we still 
lack a satisfactory understanding of the behaviour of Yang- 
Mills theories in the low energy regime. Here the coupling 
constant of the theory is large and nonperturbative effects have 
to be taken into account. 

The introduction of condensates, i.e. the (integrated) vac- 
uum expectation value of certain operators, allows one to 
parametrize certain nonperturbative effects arising from the 
infrared sector of e.g. the theory described by Q. Via the 
Operator Product Expansion (OPE) (viz. short distance ex- 
pansion), which is applicable to local operators, one can re- 
late these condensates to power corrections which give non- 
perturbative information in addition to the perturbatively cal- 
culable results. If one wants to consider the possible effects 
of condensates on physical quantities in a gauge theory, quite 
clearly only gauge invariant operators should be considered. 
The most famous example is the dimension 4 gluon conden- 
sate (^asF^y), giving rise to power corrections. Via the 
SVZ (Shifman-Vainshtein-Zakharov) sum rules |3], one can 
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extract phenomenological estimates for (<XsF^ ^ 



fjV/- 



In recent years, a great deal of interest arose in dimension 2 
condensates in gauge theories. Most attention was paid to the 
gluon condensate (A^^^ in the Landau gauge, due to the work 
of LI, as the quantity 



\ 



min/ 



mm 

uesu{N) 



^Jd^iA",f), (5) 



which is gauge invariant due to the minimization along the 
gauge orbits, could be physically relevant. In fact, as shown 
in 01 in the case of compact QED, the quantity (A^jj^) 
seems to be useful in order to detect the presence of nontrivial 
field configurations like monopoles. One can show that A^ 
can be written as an infinite series of nonlocal terms, see | 
and references therein, namely 



A^ . = 



A? (|3A" 



r 



^baMa; 



(6) 



Since the operator A^j^^ is nonlocal, it falls beyond the ap- 
plicability of the OPE annex sum rules, which refer to local 
operators. 

However, in the Landau gauge, 3;jA^, = 0, all nonlocal terms 
of expression (|6j drop out, so that A^^j^^ reduces to the lo- 
cal operator A^, hence the interest in the Landau gauge and 
its dimension two gluon condensate (A^^ A complication is 
that the explicit determination of the absolute minimum of A^ 
along its gauge orbit, and moreover of its vacuum expectation 
value, is a very delicate J^sue intimately related to the problem 
of Gribov copies i8liljiolllll[lJ,llJ. 

Nevertheless, some nontrivial results were proven concern- 
ing the operator A^. In particular, we mention its multi- 
plicative renormalizability to all orders of perturbation the- 
ory, in addition to an interesting and numerically verified re- 
lation concerning its anomalous dimension fT4llT^ . An ef- 
fective potential approach consistent with renormalizability 
and renormalization group requirements for local compos- 
ite operators (LCO) has also been worked out for this oper- 
ator, giving further evidence of a nonvanishing condensate 
(a j) ^ 0, which lowers the nonperturbative vacuum energy 
flal . The LCO method yields an effect ive gluo n mass squared 
m] - (a2) of a few hundred MeV (TUllllllJj. 

In 12011 . it was already argued that gauge (in)variant conden- 
sates could also influence gauge variant quantities such as the 
gluon propagator An OPE argument based on lattice simula- 
tions in the Landau gauge has indeed provided evidence that 
a condensate (A^^ could account for quadratic power correc- 



tions of the form 



reported in the running of the cou- 



pling co nstant as well as in the gluon propagator, see e.g. 
El 11 il El m lli. This OPE approach allows one to 
obtain an estimate for the soft part (A^)^^ originating from 
the infrared sector. The OPE can also be employed to relate 
this condensate to an effective gluon mass i2lll . 

The presence of mass parameters in the gluon propagator 
have also been advocated from the lattice perspective several 
times |27, 28, 29, 30, 31, 32], whilst effective gluon masses 
also found phenomenological use |33, 34]. 

A somewhat weak point about the operator A^j^^ is that it is 
unclear how to deal with it in gauges other than the Landau 
gauge. Till now, it seems hopeless to prove its renormaliz- 



ability out of the Landau gauge. In fact, at the classical level, 
adding Q to the Yang-Mills action is equivalent to add the 
so-called Stueckelberg action, which is known to be not renor- 
malizable (SSllsSl . We refer to |7] for details and references. 
As already mentioned, also the OPE becomes useless outside 
the Landau gauge for this particular operator 

In other gauges, there can be found other renormalizable 
local operators, which condense and give rise to a dynamical 
gluon mass. Next to the Landau gauge fT4lfl^fl7l fT ^ fls i. l3 ^ 
38] the maximal Abelian 1 39, 40], linear covariant [41 ,42,431] 
and Curci-Ferrari gauges 1 44. 45.1 have been investigated in 
the past. 

The relevant operators in these other gauges are however 
gauge variant, and as a consequence also the effective gluon 
mass. From this perspective, it is worthwhile to find out 
whether a gauge invariant framework might be found for a 
dynamical mass, and related to it for ^ power corrections. 

In order to have a starting point, we need a dimension 2 
operator that is gauge invariant. This necessarily implies a 
nonlocal operator, since gauge invariant local operators of di- 
mension 2 do not exist. We would also need a consistent cal- 
culational framework, which requires an action only contain- 
ing local terms. Therefore, we should find an operator that 
can be localized by means of a finite set of auxiliary fields, in 
such a way that the local gauge invariance is respected. As 
^min looks a bit hopeless from this viewpoint as it is a infinite 
series of nonlocal terms (|6}, we moved our attention instead 
to the nonlocal gauge invariant operator 



= — [ d^xF" 
VT J 



, -1 



ab 



(7) 



This operator caught already some attention in 3 dimensional 
gauge theories in relation to a dynamical mass generation 

In the foUowing sections, we shall show that the operator 
(0 can be localized, giving rise to a local, classically gauge 
invariant action. Afterwards, we discuss how to investigate 
the renormalizability of the action once quantized. Eventually, 
we need to introduce a slightly more general classical action in 
order to obtain a quantum action that is renormalizable to all 
orders of perturbation theory. In the case of vanishing mass, 
the equivalence of our action with usual Yang-Mills theories 
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can be shown. We shall point out the existence of a naturally 
extended version of the usual BRST symmetry. Before turn- 
ing to conclusions, we explicitly give various renormalization 
group functions, verifying the renormalizability at the practi- 
cal level. 



II. CONSTRUCTION OF THE ACTION AND ITS 
RENORMALIZABILITY AT THE QUANTUM LEVEL 

A. The action at the classical level 



B. The action at the quantum level 



In order to discuss the renormalizabilit y of JI H . we relied 
on a method introduced by Zwanziger in I47l l48ll . Instead of 
using Jl li with m coupled to the composite operators B"i^FI^^ 
and B'^^F^^, we introduce 2 suitable external sources Vpo^v and 
Vpa^iv and replace S,„ by 



B^r.F,':. 



(14) 



We can add the operator Q to the Yang-Mills action as a 
mass term via 



At the end, the sources Vappvix), Vap^vix) are required to at- 
tain their physical value, namely 



Sym + Sc 



with 



4 , ^'<v 



ab 



(8) 



(9) 



As we have discussed in 101, the action (|8} can be localized by 
introducing a pair of complex bosonic antisymmetric tensor 
fields, (Z?^y,B^y), and a pair of complex anticommuting anti- 
symmetric tensor fields, (G^v,G^v)' belonging to the adjoint 
representation, according to which 



V, 



OpfJV 



phys 



Op/jV 



phys 



--^ (5a/j8pv - 8crv6pAi) , (15) 



SO that (I14> reduces to Sm in the physical limit. 

From now on, we assume the linear covariant gauge fixing, 
implemented through 

% = J d\ (^^b"b" +b"d,Al + c"d,Dfc>') , (16) 

In (tJ, we wrote down a list of symmetries enjoyed by the 
action 



DBDBDGDGexp 



d'xG^DfD'jG';,,- 



1 
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im 



d'xB%DfD'^Bl„ 



d'x{B-B)l^F;:, 



(10) 



Therefore, we obtain a classical local action which reads 



Sym + Sbg + S,, 



where 

Sbg 



1 

4 
mi 



d X 



(B%DfD'i; 



•'''BL-(f„^D"'Dt'G'^ 



- d'xiB-B)" F" 



4 y " - V- -y^jv A^v 
which is left invariant by the gauge transformations 



(11) 



(12) 



Sym + Sbg + 5?/ , 



(17) 



i.e. in absence of the sources. Let us only mention here the 
BRST symmetry, generated by the nilpotent transformation s 
given by 



sc 



> rabc a h 

-J c c 



gfabc^bgc 



fiV 



: gf'''c''B' 



= b" , sb"^0, s^^O. 



gf cG^^+B^ 



'fjV 



(18) 



It turns out that one can introduce all the necessary exter- 
nal sources in a way consistent with the starting symmetries. 
This allows to write down several Ward identities by which 
the most general counterterm is restricted using the algebraic 
renormalization formalism l7l l4^ . After a very cumbersome 
analysis, it turns out that the action (II 1> must be modified to 



Sphys — Scl + Sgf 



(19) 



(13) with 



Scl = 



d'^x 



b nbc/~^c 
pv 



-m Xi [B^^B"^^ - Gp^G"p^) 



-"^%{B%-B'p,) 



2 , X"'^"' 
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pa"pa "^pa"^pa 



(20) 
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in order to have renormalizability to all orders of perturbation 
theory. We notice that we had to introduce a new invariant 
quartic tensor coupling X"''"', subject to the generalized Jacobi 
identity 



jrmanj^mhcd _^ ^rnhnj^amcd 



f 



•mcny^abmd ^ j^mdny^hcm q ^21) 



and the symmetry constraints 

yabcd ycdab 



bacd 



as well as two new mass couplings A,i and A-3. Without the new 
couplings, i.e. when A-i = 0, ^3 = 0, A,"''"^ = 0, theprevious 
action would not be renormalizable. We refer to fr, 50] for 
all the details. We also notice that the novel fields B^^, B^^, 
G^v and G^^ are no longer appearing at most quadratically. 
As it should be expected, the classical action Sd is still gauge 
invariant w.r.t. i\3i . 



physical operators shall not depend on the choice of the gauge 
parameter |49]. 

B. Existence of a "supersymmetry" when m = 



We define a nilpotent (anti-commuting) transformation 5j 



as 







— r"' 








- r" 




(22) 


5,. (rest) 


= 0. 





(25) 

Then one easily verifies that ( I25> generates a "supersymme- 



try" of the action since 



phys 



= 0. 



with 



5? = 0, 



(26) 



(27) 



III. FURTHER PROPERTIES OF THE ACTION 

A. Existence of a BRST symmetry with a nilpotent charge 

The BRST transformation (I18> no longer generates a sym- 
metry of the action Sphys- However, we are able to define a 
natural generalization of the usual BRST symmetry that does 
constitute an invariance of the gauge fixed action (I19> . Indeed, 
after inspection, one shall find that 



S^phys 



0, 
0, 



(23) 



Taking another look at the transformations s and s, respec- 
tively given by jl8> and J24L one recognizes that 



s = s + 8.5 , 
{5.,,?} = 0. 



(28) 



Since Sj is a nilpotent operator, it possesses its own cohomol- 
ogy, which is easily identified with polynomials in the original 
Yang-Mills fields {A'^,b" ,c" ,c"}. The auxiliary tensor fields, 

{B^v,B^v,G^v'^^v}' do not belong to the cohomology of 5^, 
because they form pairs of doublets (4^. This fact can be 
brought to use in the following subsection. 



with 



Srabcab 

2-' 



b" , Jb"=0. 



- sf"'"'c''G'' 

pv 1 ""^v aJ ^pv 



(24) 



Hence, the action Sphys is invariant with respect to a nilpo- 
tent BRST transformation J. We obtained thus a gauge field 
theory, described by the action Sphys, il9i . containing a mass 
term, and which has the property of being renormalizable, 
while nevertheless a nilpotent BRST transformation express- 
ing the gauge invariance after gauge fixing exists simultane- 
ously. It is clear that stands for the usual BRST transforma- 
tion, well known from literature, on the original Yang-Mills 
fields, whereas the gauge fixing part Sgf given in ilOi can 
be written as a 7-variation, ensuring that the gauge invariant 



C. Equivalence with Yang-Mills gauge theory when m = 

If the mass m = 0, we would expect that the action ( I19l l 
would be equivalent with the usual Yang-Mills gauge the- 
ory, since in the nonlocal formulation (|9jl, we would have in- 
troduced "nothing". In the local renormalizable formulation 
( II9I 1. this would also be trivially true when X"'"^'^ = as then 
we would only have added a -although quite complicated- 
unity to the Yang-Mills action. Unfortunately, since renor- 
malization forbids setting X"^"^ = 0, we must find another ar- 
gument to relate 5'™,=^ to the usual Yang-Mills gauge theory. 
As proven in 150.] . the "supersymmetry" 5s of j25> can be used 
to show that 



{G„{xi ,.. . ,x„)) 



SyM+Sgf 



— {Gn{xi , • • • ,Xfi))„m= 



phys 



(29) 



where 



Gn{x\,...,Xn) = A(x\) . . .A{xi)c(xi+i) . . .c(xj)c{xj+\) . . .c(xk)b(xk+\) . . .b(xn) 



(30) 
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is a generic Yang-Mills functional. The expectation value of 
any Yang-Mills Green function, constructed from the fields 
{A^jC'^^Tfjb"} and calculated with the original (gauge fixed) 
Yang-Mills action Sym + Sgf, is thus identical to the one cal- 
culated with the massless action S' 



1=0 
phys ' 



where it is of course 
assumed that the gauge freedom of both actions has been fixed 
by an identical gauge fixing. 

The foregoing result also reflects on the renormalization 
group functions. As usual, we employ a massless renormal- 
ization scheme known as the MS scheme. As a consequence, 
we can set m = in order to extract the ultraviolet behaviour 
Using ( I29> , we conclude that all the renormalization group 
functions of the original Yang-Mills quantities are not affected 
by the presence of the extra fields or couplings. This fact shall 
be explicitly verified in the next section. 



IV. EXPLICIT RENORMALIZATION AT TWO LOOP 
ORDER 

Having proven the renormalizability of the action il9\ . we 
shall now compute explicitly the two loop anomalous dimen- 
sion of the fields and the one loop p-function of the tensor 
coupling A,"''"'. The corresponding details can be found in ff\ 
for one loop results, while two loop results are discussed in 

We have regarded the mass operator as an insertion and split 
the Lagrangian into a free piece involving massless fields with 
the remainder being transported to the interaction Lagrangian. 
To renormalize the operator, we insert it into a massless Green 
function, after the fields and couplings have been renormal- 
ized in the massless Lagrangian. An attractive feature of the 
massless field approach is that we can use the MINCER al- 
gorithm to perform the actual computations. This algorithm, 
15111 . written in the symbolic manipulation language Form, 



is devised to extract the divergences from massless 
2-point functions. The propagators of the massless fields in an 
arbitrary linear covariant gauge are 



{c"{p)c''{-p)) 



{Kvip)BU-p)) 



{CMgU-p)) = 



5"* 

1t2 



8,v - il-a)^ 



- -j^ [§^o8vp - 8/jp8vo] , (31) 



where p is the momentum. The necessary Feynman diagrams 
were generated automatically with Qgraf 15411 . 

We first checked that the same two loop anomalous dimen- 
sions emerge for the gluon, Faddeev-Popov ghost and quarks 
in an arbitrary linear covariant gauge as when the extra local- 
izing fields are absent. It was also explicitly verified that the 
correct coupling constant renormalization constant is found. 
These results are in agreement with the general argument of 
the previous subsection. 

We have implemented the properties (12 1> and \22\ of the 
y^abcd coupling in a Form module, while it was assumed that 



1' 



acde 1 bcde 



1 



~ Na 



(32) 



which follows from the fact that there is only one rank 2 in- 
variant tensor in a classical Lie group. 

At two loops in the MS scheme, we find that 



yB(a,A.) = YG(aA) = (a-3)fl 



^ 61 

h2a 

4 6 
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cbd 



(33) 



r 



where Na is the dimension of the adjoint representation of the 

2 

colour group, a = -jj^ and we have also absorbed a factor of 
^ into X"''"' here and in later anomalous dimensions. These 
anomalous dimensions are consistent with the general obser- 
vation that these fields must have the same renormalization 
constants, in agreement with the output of the Ward identities 
Q]. A check on ( I33> is that after the renormalization of the 
3-point gluon B"^y vertex, the correct gauge parameter inde- 
pendent coupling constant renormalization constant emerges. 

We also determined the one loop p-function for the X"'"^'' 
couplings. As this is present in a quartic interaction it means 



that to deduce its renormalization constant, we need to con- 
sider a 4-point function. However, in such a situation the 
Mincer algorithm is not applicable since two external mo- 
menta have to be nullified and this will lead to spurious in- 
frared infinities which could potentially corrupt the renormal- 
ization constant. Therefore, for this renormalization only, we 
have resorted to using a temporary mass regularization intro- 
duced into the computation using the algorithm of |55] and 
implemented in Form. Consequently, we find the gauge pa- 
rameter independent anomalous dimension 
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I 



(34) 



from both the X"'^"'B);^B^''^B^pB''°P and 

V^"'Bl^B^P''G'^pG''''f' vertices where df"^ is the totally 
symmetric rank four tensor defined by 



"A 



bed 



Tr 



/ rparpi^ rpC rpd 

yA^A ^A^A 



(35) 



with denoting the group generator in the adjoint represen- 
tation, 1 56]. Producing the same expression for both these 
4-point functions, aside from the gauge independence, is a 
strong check on their correctness as well as the correct im- 
plementation of the group theory. Moreover, as it should be, 
^abcd gjjjoys the same symmetry properties as the tensor A,"'"^^, 
summarized in (I22> . 

We notice that A,"^'^"^ = is not a fixed point due to the extra 



J 



A-"''"' -independent terms. If we had not included the X"*^''- 
interaction term in the original action, then such a term would 
inevitably be generated at one loop through quantum correc- 
tions, meaning that in this case there would have been a break- 
down of renormalizability. 

Finally, we turn to the two loop renormalization of the mass 
m. The corresponding operator can be read off from M9\ and 
is given by 



M = {Bl. 



pv 



R )F" 

"pV J ^ pv ■ 



(36) 



We insert 5V/ into a A'^-B^^ 



2-point function and deduce 
the appropriate renormalization constant, leading to the MS 
anomalous dimension 



yo(a,X) = 



2\-TpNf- 



a - (iTpNfCA+^TpNfCF - ^C; 



ahe rcdei adhc ^ 



r^"'r'"t37) 



IIMa 



r 



as the two loop MS anomalous dimension. The gauge param- 
eter independence is again a good check, as the operator 9vt is 
gauge invariant. 



V. CONCLUSIONS 

We added a nonlocal mass term (|9jl to the Yang-Mills ac- 
tion Q, and starting from this, we succeeded in construct- 
ing a renormalizable massive gauge model, which is gauge 
invariant at the classical level and when quantized it enjoys 
a nilpotent BRST symmetry. This BRST symmetry ensures 
that the expectation value of gauge invariant operator is gauge 
parameter independent. We have also proven the equivalence 
of the massless version of our model with Yang-Mills gauge 
theories making use of a "supersymmetry" existing between 
the extra fields in that case. We presented explicit two loop 
renormalization functions, thereby verifying that the anoma- 
lous dimensions of the original Yang-Mills quantities remain 
unchanged. 

Many things could be investigated in the future concerning 
the gauge model described by M9\ . 

At the perturbative level, it could be investigated which 
(asymptotic) states belong to a physical subspace of the 
model, and in addition one should find out whether this phys- 
ical subspace can be endowed with a positive norm, which 



would imply unitarity. The nilpotent BRST symmetry MAI 
might be useful for this. 

The model M9\ is also asymptotically free, implying that 
at low energies nonperturbative effects, such as confinement, 
could set in. Proving and understanding the possible confine- 
ment mechanism in this model is probably as difficult as for 
usual Yang-Mills gauge theories. 

It would also be interesting to find out whether a dynami- 
cally generated term m{B — B)F might emerge, which in turn 
could influence the gluon Green functions. This might also 
be relevant in the context of gauge invariant power correc- 
tions, an issue that recently has also attracted attention from 
the gauge/string duaUty side, the so-called AdS/QCD iSTllssil . 
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